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ABSTRACT 

The X-ray background can be used to constrain the accretion history of Supermassive Black 
Holes (SMBHs) in Active Galactic Nuclei (AGN), with the SMBH mass density related to 
the energy density due to accretion. A knowledge of the hard X-ray bolometric correction, 
K2-iokoV is a vital input into these studies, as it allows us to constrain the parameters of the 
accretion responsible for SMBH growth. Earlier studies assumed a constant bolometric cor- 
rection for all AGN, and more recent work has suggested accounting for a dependence on 
AGN luminosity. Until recently, the variations in the disk emission in the UV have not been 
taken into account in this calculation; we show that such variations are important by construc- 
tion of optical - to - X-ray SEDs for 54 AGN. In particular, we use FUSE UV and X-ray 
data from the Uterature to constrain the disk emission as well as possible. We find evidence 
for very significant spread in the bolometric corrections, with no simple dependence on lu- 
minosity being evident. Populations of AGN such as Narrow-Line Seyfert 1 nuclei (NLSls), 
Radio Loud and X-ray Weak AGN may have bolometric corrections which differ systemati- 
cally from the rest of the AGN population. We identify other sources of uncertainty including 
intrinsic extinction in the optical-UV, X-ray and UV variability and uncertainties in SMBH 
mass estimates. Our results suggest a more well-defined relationship between bolometric cor- 
rection and Eddington ratio in AGN, with a transitional region at an Eddington ratio of ^ 0.1, 
below which the bolometric correction is typically 15 — 25, and above which it is typically 
40 — 70. We consider the potential implied parallels with the low/hard and high/soft states in 
Galactic Black Hole (GBH) accretion, and present bolometric corrections for the GBH binary 
GX 339-4 for comparison. Our findings reinforce previous studies proposing a multi-state 
description of AGN accretion analogous to that for GBH binaries. Future calculations of the 
SMBH mass density may need to take into account the possible dependence of K2-iokcV on 
Eddington ratio. 
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1 INTRODUCTION 

Numerous authors have investigated the relationship between past 
Active Galactic Nuclei (AGN) activity and the local supermassive 
black hole (SMBH) population. A detailed study of the relation- 
ship between the two can yield important insights into the nature 
of AGN activity, providing us with parameters such as the typical 
accretion efficiency and Eddington ratio for AGN. 

Early work by Soltan C198Z) demonstrated how the total mass 
density of SMBHs in the local universe could be inferred from 
counts of quasars/ AGN (hereafter refer red to as the 'Soltan ar- 
gument'). Shankar et al. (2004) and Marc oni et al. I f2004) demon- 
strate how the mass function of local SMBHs can be matched with 



that expected from accreting AGN, for appropriately chosen accre- 
tion efficiencies and Eddington ratios. 

Various observational studies have convincingly established 
that the X-ray Background (XRB) can be resolved as the sum of 
the spectra from AGN (e.g. as summarised in Hasinger 2004). It 
is then possible to use the X-ray background light as a 'signature' 
for accretion causing SMBH growth, and so the energy density in 
the XRB can be converted into a mass density of SMBHs expected 
in the present epoch, taking into account the efficiency of SMBH 
ac cretion in the p rocess. The basic method for this is summarised 
bv lFabiadfcooai) : if the observed energy density due to accretion 
is given by and the mass-to-radiation conversion efficiency is 
77, then the local SMBH mass density pen is: 
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where (z) is the mean redshift of accreting AGN. 

However, we require a bolometric correction K2-iokcV to 
scale up the 2-10 keV XRB spectral intensity to obtain the total 
accretion energy density: 



4tt 



(2) 



where /xrb is the 2-10 keV XRB spectral intensity. The quantity 
Ut can then be used in equation l[T) to determine the SMBH mass 
density. Despite previous detailed studies into AGN and quasar 
spectral en e rgy di stributions (SEDs) such as the seminal work by 
lElvis et alj ( 1 1994) . there remain significant issues regarding the 
determination of this bolometric correction. This can have wide- 
ranging effects, since the choice of bolometric correction will di- 
rectly constrain the allowed efficiencies for which the Soltan argu- 
ment can be satisfied. The efficiency of accretion is directly related 
to black hole spin, and so the quantity K2-iokcV needs to be well 
understood. 

Previously, bolometric corrections have been determined from 
a Mean Energy Distributio n (MED) for AGN calculated from 47 
quasars telvis et al] Il994h . It is important to note that the in- 
frared (IR) emission is known to be re-processed from the ultra- 
violet (UV), as discussed by Antonucci ( 1993). The bolometric 
corrections of Elvi s et al.l | |1994|) . which include the IR emission 
in the bolometric luminosity, effectively count part of the emission 
twice and the resulti ng bolometric corrections would be too high 
jMarcoiri et alj[2004h . and would overestimate the accretion lumi- 
nosity. Another key issue is the configuration of the absorbing ma- 
terial responsible for producing the IR emission: while the UV and 
X-ray emission from the disk and corona are thought to be emitted 
more isotropically, the IR emitting material is unlikely to cover the 
source completely. Therefore in the process of converting the IR 
flux into a luminosity, the implicit assumption of isotropic emis- 
sion would yield an artificially high IR contribution. To recover 
the intrinsic accretion luminosity in AGN (which emerges in the 
optical, UV and X-ray bands), in general the IR emission should 
be excluded when calculating bolometric corrections. There could 
be an important sync hrotron contribution t o the IR in Radio Loud 
AGN as discussed in IMarconi et J] j2004h . but we do not explore 
this further here (on ly around 10 per cent of our objects are radio 
loud). Additionally, felvis et alj ( IT994h themselves caution against 
using the MED to determine global properties of quasars and AGN, 
since the individual SEDs used to determine the MED themselves 
show significant variation over many wavelength bands. 

More recent studies have accounted for variations in bolo- 
metric correction using a re-scalable template SED. Mar coni et al.l 
( l2004h account for variations in AGN SEDs by using the well 
known anti-correlation between t he optical-to-X-ray spectral in- 
dex (aox) and 2500 A luminosity l lVignali et alj2003h to construct 
such a template, where aox is defined as: 



aox 



log[L42500A)/L42keV)] 



(3) 



log[i/(2500A)///(2keV)] 

IMarconi et aP ( |2004|) renormalise their template SED to a 
particular value of aox to obtain the variation of bolometric 
correction with lum inosity. A similar approach is followed by 
iHopkins et al.l j2006h . but employs a template SED generated from 
the averages of real SEDs in different wavebands. The variations of 
bolometric correction with luminosity calculated by both of these 
studies are reproduced in Fig. [T] The latter paper does advise that 
intrinsic spread in AGN SEDs could give rise to variation of a factor 




10* 10" 10" 10-' 10-" 10-" lO" 
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Figure 1. Bolometric correction as a function of luminosity, with spread 
calculated from Mon te-Carlo realisations of a tem plate AGN SED, from 
IMarconi et"ZI j2004l) (lower grey filled area) and iHopkins et"ai] j2006l) 
(higher red filled area). 



of ~ 2 in the bolometric correction, even when luminosity depen- 
dence is taken into account. 

We have constructed basic SEDs for 54 AGN and calculated 
their bolometric corrections, for comparison with the trends sug- 
gested by these authors. Historically it has been established that 
the optical to X-ray SED of AGN can be relatively well described 
by a black body disk spectrum peaking in the UV and a power- 
law model extending to hard X-rays; but most emission from lower 
luminosity objects is unobserved by virtue of it peaking in the ab- 
sorbed far ultra-violet (FUV). We therefore used UV observations 
as the foundation for our study, to constrain the disk emission com- 
ponent. We employed data from a range of sources for optical and 
X-ray SED points and fitted all the available multi-wavelength data 
to a basic disk and power-law model. 

The comparison between accretion on large and small scales 
is a field of much interest, and recent work has sought to describe 
accretion across all mass scales in terms of similar p hysical pro- 
cesses and accretion states (e.g. iMcHardv et al]|2006l) . The accre- 
tion states in Galactic Black Hole binaries (GBHs) are seen to be 
separated by different Eddington ratios (Lboi/LEdd), and to facil- 
itate comparison with AGN we present the bolometric corrections 
as a function of Eddington ratio. We compare our findings qual- 
itatively to the current understanding of state transitions in GBH 
X-ray binaries. 



2 DATA SOURCES 

The starting point of our study is the sample of 85 low redshift 
(z < 0.7) AGN surveyed by the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE) as presented bv lScott et alj 12004 ): this provides a 
measure of the "big blue bump" (BBB) emission from the ther- 
mal accretion disk in the UV. Observations from HST-STIS have 
been used where FUSE data were not available. Optical and X- 
ray data have been gathered from publications in the literature em- 
ploying a variety of instruments, including the Hubble Space Tele- 
scope (HST), Kitt Peak National Observatory (KPNO), the Ad- 
vanced Satellite for Cosmology and Astrophysics (ASCA, as pre- 
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sented in the TARTARUS Databas^B) and XMM-Newton. Complete 
datasets (UV, X-ray and optical if available) were obtained for 54 
AGN of which 48 were drawn from the FUSE sample of Scott et 
al. The data sources used are listed in Table [All in Appendix lAl 



2.1 Sample Selection 

The study o f lElvisetal.lJl994h is consists primarily of X-ray bright 
AGN drawn from an X-ray selected sample.Therefore an X-ray 
bolometric correction calculated from their mean SED would be 
an underestimate for the AGN with more canonical X-ray lumi- 
nositites. Si nce our sample d raws largely from the AGN sample 
presented in IScott et al.U2004h . any selection biases present in their 
sample will necessarily transfer directly to our study. By nature, the 
FUSE study will be UV bright and indeed contains few very low 
luminosity (Lx < lO^^ergs^^) Seyferts, and so for the lowest lu- 
minosity objects in table|A| we obtain UV data from other sources. 
The preponderance for brighter Seyferts and quasars should be 

noted at the outset. 

Prior to the work of'Scott et al.l j2004h . lTelfer et all ( l2002h pro- 
vided an analysis of the UV spectral properties of 184 QSOs with 
UV observations from HST. They found a distinctly softer average 
Extreme-UV (EUV) spectral slope (qeuv = —1.76 ± 0.12 com- 
pared to Scott et al.'s oeuv = -0.56ln for ,U oc v"" '"^). 
Their sample is brighter in the UV than the lScott et al.l ( l2004h sam- 



ple, with an ave rage valu e of log(AL^^^^^) ~ 45.9 compared to 
~ 45.0 for the IScott et a l. (2004) sample. Additionally, the me- 
dian redshift of the Telfer et al., t2002) sample is 0.96 compared 



to 0.1 for the lScott et alj j2004h sample. By cosmic downsizing, it 
is reasonable to expect significantly higher SMBH masses for the 
higher redshift Telfer et al. sample. This systematic difference in 
the average mass could account for the difference in spectral shape, 
possibly reducing the fraction of luminosity from the BBB. If the 
SEDs from the Telfer et al. sample were assumed to lie on the aox 
- 2500A luminosity relation, we would expect similar bolomet- 
ric corrections, with the possibility of some reduction due to the 
smaller disk fraction extrapolated from th e softer EUV slope . 

Without a record of the objects in the lTelfer et al. I Jiooih sam- 
ple along with mass estimates (to constrain the shape of the disk 
model in the UV), it is difficult to make any predictions as to how 
the bolometric corrections might vary between these two samples. 
We concentrate on the lower redshift objects of interest here, as 
presented in the FUSE sample. 

When attempting to put together this sample, the prime con- 
cern was to assemble a catalogue of objects for which at least UV, 
X-ray and mass data were available and so the sample selection has 
consequently been constrained by the availability of published data 
on the sources of interest. However the subsequent analysis reveals 
that despite the biases of the parent samples, there are a range of 
object classes represented, notably both radio loud and radio quiet 
AGN as well as a substantial number of Narrow Line Seyfert 1 nu- 
clei (NLS Is). 



3 CONSTRUCTION OF SEDS 

The multi-wavelength data points were converted into XSPEC 
spectra files using the FLX2XSP utility, part of the FTOOLS pack- 
age. The data points (corrected for redshift) were fit to the DISKPN 
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model in the XPSEC package, with a simple power law at X-ray 
energies up to 250keV (POWERLAW model) or, if the X-ray pho- 
ton index was steeply decreasing, a broken power law increasing 
steeply within the UV disk region (BKNPOWER model) in order 
to prevent the extrapolation of the powerlaw from producing arti- 
ficially high fluxes at low energies. The normalisation of the disk 
component was constrained using SMBH mass estimates available 
in the literature and was kept frozen during fitting. The disk tem- 
perature was seeded with a typical value of O.OlkeV and subse- 
quently allowed to vary, and the inner radius of the accretion disk 
was frozen at 6GM/c^; appropriate for the high accretion efficien- 
cies found from previous applications of the Soltan argument. 

If multiple epoch X-ray observations were available, the mean 
2 — lOkeV flux and photon index from the available observations 
were used. The TARTARUS database provides both an 'absorbed' 
and an 'unabsorbed' model fit to the ASCA X-ray data, and if the 
fit parameters (F2-iokcV, 2-10 keV flux and photon index F) were 
significantly different between the two fits, the fit parameters with 
lower reduced chi-squared were used. If both had similar reduced 
chi-squared values but different model fit parameters, the phys- 
ically simpler 'unabsorbed' parameters were used. Where errors 
were not provided (e.g. on individual fluxes and photon indices in 
the Tartarus database), errors of ±0.1 on spectral indices and 
±10 per cent on fluxes w ere assumed. Optical points were taken 
primarily from the study of lBaskin & LaoiH2005h and NAS A/IPAC 
Extragalactic Database (NED) photometry. The optical points did 
not always follow the shape implied by a BBB expected in the 
optical-UV. When the optical fluxes lay significantly above the UV 
points, deviating markedly from the best disk model fit, they were 
not included in a revised SED fit. This optical excess may be due 
to host galaxy starlight contamination, or significant reddening in 
the UV causing a real decrease in the UV fluxes with respect to the 
optical fluxes. In this study we make the assumption that the AGN 
have at least a basic BBB shape, in order to fit a disk model to the 
UV-optical region. 

The hard X-ray bolometric correction for each AGN was 
determined by dividing the luminosity in the energy range 
0.001 — 250keV (i.e. from optical wavelengths to hard X-rays) by 
the 2-10 keV luminosity. Some sample SEDs generated via this 
process are shown in Fig. |2l with bolometric correction plotted 
against luminosity for the 54 AGN in Fig. [3] The bolometric cor- 
rections for the AGN in the sample are reproduced in Table [T] 



3.1 Checking SED shapes against the Literature 

The anti-correlation between qqx and 2500A luminosity dis- 
cussed in Section [T] usefully provides a constraint on the physical 
processes at work, as it tells us that the fraction of power in the 
accretion disc corona (e mitted in the X-rays ) decreases with in- 
creasing UV luminosity. IStrateva et alj ( l2005h used 228 optically 
selected AGN across a range of redshifts (0.01 < z < 6.3) to pro- 
vide a robust characterisation of this relation. They confirm that the 
Qox — L„ (2500A) relation represents an underlying physi cal cor- 
relatio n rather than an accumulation of systematics. Steff en et al.l 
( l2006i) extend their work to a larger range in luminosities by includ- 
ing fainter AGN (from the Extended Chandra Deep Field-South) 
and more luminous AGN (from the Bright Quasar Survey). This 
provides the most recent and robust determination of the relation, 
spanning the largest space in the luminosity-redshift plane. We 
comp are the values of q qx generated from our SEDs to the best 
fit of lSteffen et alj ( |2006|) in Fig.g] 

It can be seen that the values for aox obtained lie within the 
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Table 1. Bolometric corrections for all AGN in the sample. [1] 2-10 keV luminosity from XSPEC fit. May differ slightly from the values in the source literature 
due to slight variations in the best fit parameters. [2] Bolometric accretion luminosity is the total luminosity in the energy range 0.001 — 250keV. [3] 2-10 
keV bolometric correction. Errors on bolometric corrections are 'max/min'; i.e. they represent the range of values achievable when the input parameters are 
varied within their error limits. [4] Logarithm of Black Hole mass Mbh in units of solar mass. *The mass estimate for Mrk 506 was taken from a study with 

systematically higher masses than in other studies used here, so the mass estimate for Mrk 506 has been scaled down appropriately. [5] "NLSI" denotes a 
Narrow Line Seyfert 1 nucleus, "RL" denotes a radio loud AGN and "XRW" denotes an X-ray weak AGN. 
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Figure 3. Hard X-ray Bolometric Con'e ction ag ainst 2-10 keV luminosity for all AGN in sample. The grey (lower) and red (higher) shaded areas represent the 
bolometric corrections (with spread) of lMarconi et al. (2004) and Hopkins et al. (2006) respectively. The dashed black and red lines represent the bolometric 
corrections that would have been obtained by [Marconi et alj ( i2004.) and .Hopkins et ah t2006) if the IR had been included at a level of ^ 1/3 of the total 
luminosity. Blue empty squares represent Narrow Line Seyfert 1 nuclei (NLSls), black empty circles represent Radio Loud AGN (RL) and red diamonds 
represent X-ray weak. The AGN in the sample: 1 - Ark 564, 2 - I Zw 1, 3 - PG 1116-H215, 4 - KUG 1031-H398, 5 - PG 1512-1-370, 6 - Mrk 335, 7 - PG 
1322-1-659, 8 - TON 951, 9 - PG 0026-H129, 10 - PG 1402-1-261, II - Mrk 478, 12 - PG 1211-1-143, 13 - Ark 120, 14 - TONS 180, 15 - 4C -H34.47, 16 - PG 
0052-1-251, 17 - PG 0804-H761, 18 - PG0947-^396, 19 - PG 0953-H414, 20 - PG 1114-1-445, 21 - PG 1115-H407, 22 - PG 1216-H069, 23 - PG 1307-1-085, 24 - PG 
1415-1-451, 25 - PG 1444-1-407, 26 - PG 1626-1-554, 27 - PG 1229-H204, 28 - Mrk 205, 29 - Mrk 79, 30 - Mrk 279, 31 - MR 2251-178, 32 - 3C 273, 33 - Mrk 
1383, 34 - Fairall 9, 35 - NGC 985, 36 - 3C 382, 37 - PG 2349-014, 38 - NGC 4395, 39 - NGC 4593, 40 - NGC 3516, 41 - NGC 3227*, 42 - NGC 7469, 43 - 
NGC 6814, 44 - PG 1011-040, 45 - NGC 3873, 46 - Mrk 290, 47 - PG llOO-t-772, 48 - PG 1352-1-183, 49 - ESO 141-G55, 50 - PG 1048-H342, 51 - IIZw 136, 
52 - Mrk 509, 53 - Mrk 506, 54 - Mrk 876. * De-reddened UV fluxes assumed for NGC 3227 to allow a fit to flie disk model - see lCrenshaw et al.l i200ll) for 
discussion. 



limits of the lSteffen etaP (l2006h fit, implying sensible reconstruc- 
tion of the SED shapes. 



4 DISCUSSION OF SYSTEMATICS 

Fig. [5] shows that the bolometric correct ions obtained do not con- 
vincingly follo w the forms suggested by ^ Marconi et al. and 
iHopkins et al.l ([2006); the reasons for this spread could be diverse, 
ranging from selection effects to intrinsic properties of the AGN 
SED. In this section we outline a number of systematic effects 
which may affect the bolometric correction for individual AGN. 



4.1 Intrinsic Exctinction 

The shape of the BBB in the optical-UV region of AGN spec- 
tra is notoriously prone to effects of reddening. Without properly 
correcting for de-reddening effects, it is impossible to know the 
intrinsic emission of AGN in these wavebands and therefore it 
becomes difficult to constrain energy output from AGN. Galac- 
tic extinction has been investigated by numerous authors and val- 
ues of the Galactic reddening in the standard form E{B — V) are 
readily available for the AGN in our sample, along with standard 
Galactic extinctio n cur ves. All the observ ations presented in the 
IScott et all ( l2004h and IShang et al.l l l2005h samples have already 
been de-reddened for Galactic extinction according to the extinc- 
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Figure 2. Some sample SEDs. Filled circles are the data points obtained 
from the sources detailed in section|2]and solid lines show the model fits. 



tion curve of lCardelli etal] (Il989h with values for Galactic extinc- 
tion provided bv lSchlegel et alj ( ll998h . 

However the process of accounting for intrinsic extinction due 
to the AGN itself is more problematic, and such reddening can 
vary significantly between different AGN. There remains much 
debate over the shape of the wavelength dependence of the ex- 
tinction. Forms which have been proposed included Small Magel- 
lanic Cloud-type (SMC) extinction, Milky Way-type exction (both 
derived by Pe i 1992) , or a curve which is flatter in the far UV 
( iGaskell et al1[2bo4h . [Hopkins et alj ( 120041 find that an SMC-type 
extinction curve best accounts for the reddening in a sample of 
Sloan Digital Sky Survey (SPSS) quasars, in preference to Mikly 
Way-type or iGaskell et al ] ( l2004l) -type reddening. They identify 
that ~ 91 per cent of the sample has reddening E{B ~V) < 0.055 
(taking into account the effects of quasars being 'extincted out' 
due to unusually high reddenings). To obtain a relatively con- 
servative estimate of the effects of SMC-type reddening on the 
bolometric correction, we apply this reddening curve along with 
E{B — V) = 0.055 to th e 12 AGN from our sample coincident 
with the lShang et alj ( l2005h (for which the optical~UV coverage is 
more complete than the rest of our objects). We find that applying 
such de-reddening produces an increase in «:2-iokcV of a factor of 
~ 2.0. A sample of the SEDs generated including de-reddening 
corrections fo SMC-type extinct ion are reproduced in Fig . |5] 

In as yet un published work.lOaskell & Benken ( 120051) used the 
17 AGN in the Ishang et al] j2005h sample to produce a revised 
mean extinction curve for AGN, by comparing each AGN SED 
with three template 'low-reddening' objects: 3C 273, PG 1 100-1-772 
and PG 0953-1-414. These templates cover three different classes of 
AGN (namely radio loud, lobe dominated and radio quiet) thereby 
hoping to average out at least some of the effects of intrinsic object- 
to-object spectral variations in determining the average optical-UV 
intrinsic extinction curve. In the process of normalising their ex- 
tinction curves, they obtain E{B — V) values for each of the AGN 
in their sample. 
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Figure 4. Plot of aox against 2500A monochromatic luminosity for the 
AGN in our samp le. Key as in Fig. [3] The solid line represents the best fit to 
the Steffen et al. (2006) sample, with the dashed line representing the best 
fit from the study of Strateva et al. (2005). The dotted line represents the 
best fit from Vienali et al. (2003), provided for reference and comparison 
with Fig. 1 1 of Strateva et al., dlOOS^ . Grey shaded area gives the spread in 
the Steffen et al. best fit. 
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Figure 5. Comparison of a sample of SEDs with and without de -reddenin; 
due to intrinsic extinction (using the SMC extinction curve of |Pei[l99^ ^ 
assuming a reddening of i5(i? — y) = 0.055. Filled black circles and black 
solid lines are the data and model fit to the pre-dereddened data respectively, 
whereas empty blue squares and blue dashed lines are the data and model 
SED fits to the de-reddened data points. 
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Figure 6. Comparison of a sample of SEDs with and wi thout de-reddening 
due to intrinsic extinction (using the extinction curve o f iGaskell & Beniced 
l2005h . Key as in Fig.|5] 



We also used their mean extinction curve along with the values 
ofE{B- V) specific to each AGN from Gaskell & Benker ( 20q1) 
to construct de- r edden ed SEDs for the 12 AGN coincident with 
the Ishang et al ] (l2005h sample. A sample of these are plotted in 
Fig. |6l with the resulting bolometric corrections given in Fig. [7] 
An increase in bolometric correction from de-reddened SEDs is 
again clearly visible; the average increase is a factor of ^ 1.9, very 
similar to that for SMC reddening. The use of individual values of 
E[B — V) for each AGN also hints at an increase in reddening 
at lower luminosities. In summary, we estimate that the maximum 
effect of intrinsic reddening is to increase «:2-iokcV by a factor 
~ 2. 

An overarching study of reddening mechanisms in AGN, cov- 
ering the luminosity range from Quasars to Seyferts would be ex- 
tremely illuminating and allow considerable progress to be made 
in understanding the AGN SED. Nonetheless, it is clear from our 
results that the effect of optical-UV reddening will be to increase 
the X-ray bolometric correction, which will have a well defined 
effect on the SMBH mass density estimates from the XRB. If de- 
reddening uniformly increases bolometric corrections as seen here, 
then we require a greater average mass-to-radiation efficiency, r) to 
satisfy the Soltan argument. This would imply more rapidly spin- 
ning black holes and reduces the relative importance of radiatively 
inefficient modes of accretion in the X-ray background. 



4.2 Narrow Line Seyfert 1 Nuclei, Radio Loud AGN and 
X-ray Weak AGN 

Narrow Line Seyfert 1 (NLSl) nuclei are well known to exhibit 
spectral complexity at X-ray energies; this manif ests as sharp spec- 
tral drops or gradual curvature in the spectrum l lGalldl2006l) . The 
common paradigm for NLSl nuclei assumes a central black hole 
of lower mass than that typical for other classes of AGN, with ei- 
ther ionized reflection or partial covering givin g rise t o the char- 
acteristic narrow lines in the optical spectrum. iGalld ( l2006h also 
discusses the likelihood that NLSls could represent the nucleus in 
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Figure 7. Comparison of bolometric corrections with and without de- 
reddening due to intri nsic extinction (using the extinction curve of 
iGaskell & Benkeill2005l) . The filled black circles represent bolometric cor- 
rections before de-reddening, the open blue squares give the corresponding 
corrections after de-reddening. Other conventions as in Fig. |3] The AGN 
in this subsample: 1 - PG 1322+659, 2 - TON 951, 3 - 4C +34.47, 4 - PG 
0052+251, 5 - PG 0947+396, 6 - PG 0953+414. 7 - 3C 273, 8 - PG 2349- 
014, 9 - Mrk 290, 10 - PG 1 100+772, 1 1 - Mrk 509, 12 - Mrk 506. 



either a low or normal flux state, and that NLS 1 nuclei could switch 
between the two on typical X-ray variability timescales for AGN. 
We note in Fig. [3] that a large number of bolometric corrections 
for NLS 1 nuclei are seen to lie systematically above the form sug- 
gested byl Marconi et al.l ( l2004h . NLSl nuclei seem to represent an 
interesting class of object which have a larger spread in bolomet- 
ric correction than 'normal' AGN, with a possible marginal skew 
towards higher bolometric corrections. 

We have hitherto not made any consideration for the radio 
properties of these AGN, despi te the known links between their 
X-ray and radio emission (Shast riet all 19931) . Strongly radio loud 
objects potentially present complications when determining the X- 
ray emission d h'ectly associated with accretion, as pointed out by 
IStrateva et al.l ( l2005h who exclude this class of objects from the 
determination of the qqx — Ly(2500A) relation. This complica- 
tion arises because the observed X-ray flux may contain com- 
bined emission from both accretion and jet processes ( Dermer et al.l 
fT997h ; since we wish to calculate intrinsic accretion luminosities 
using the bolometric corrections, it is necessary to account for the 
energy output of the jet or exclude these objects altogether. Ac- 
counting for the X-ray flux due to jet emission is problematic, as 
it is unclear what fraction of X-ray emission in radio loud AGN is 
from jets. Fortunately, we are able to gain an estimate of the impor- 
tance of this effect with one of our radio loud AGN, 3C 273 thanks 
to the recent publication of XMM observations taken at a historic 
jet minimum by Tiirler et al. ( 2006). Using their 2-10 keV flux and 
photon index instead of ASCA data, we obtain a bolometric correc- 



tion of 2QAZ2X) compared to the value l-i.lZta obtained using the 
ASCA flux. A larger bolometric correction is in better concordance 
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with the expected value from iMarconi et al] | |2004|) . but still falls 
significantly short (the expected bolometric correction at this lumi- 
nosity is between 150-200). Thi s may imply some significant jet 
contribution still remaining in the lTiirler et al j20 0^ flux, or the low 
bolometric correction for 3C 273 may simply be another example 
of the large spread we see in the rest of our results. 

The curious object PG 1011-040 has a bolometric correction 
of ~ 10 00, the highest observe d in our sample. This AGN is iden- 
tified bv' Gallagheretal.1 ( 1200 ih as a member of the 'X-ray weak' 
AGN class, which exhibit particularly low ratios of soft X-ray flux 
to optical flux. Our inclusion of multiwavelength data confirms this 
with the high bolometric correction. These authors discuss the pos- 
sibility of these unusually low X-ray fluxes being due to either in- 
trisic X-ray weakness or a high level of X-ray absorption. They 
note that PG 1011-040 displays no spectral evidence of significant 
X-ray absorption or variability, implying that the high bolometric 
correction we obtain is probably typical for the object and does not 
represent a snapshot of an unusual state. They also point out that its 
low H/3 linewidth would place it in the NLSl subclass. 

In general, a si gnificant proportion of AGN may be X-ray 
weak, and a study bv lBrandt et al.l KOOOf) identifies ^11 per cen t 
of the optically selected QSO sample of i Boroson & Green j 19921) 
(the Bright Quasar Survey) as falling into this category. If a simi- 
lar figure holds for the general AGN population, the presence of a 
substantial proportion of AGN with bolometric corrections of order 
^ 1000 could significantly affect the SMBH density calculated via 
the Soltan argument. 

Radio loud and X-ray weak AGN are known to exhibit devia- 
tions from the qqx — L„(2500A) relation. NLSl nuclei may have 
unusual bolometric corrections as a result of their spectral com- 
plexity. It is evident from Fig. [3] that these objects also represent 
many of those which deviate significantly from the expected depen- 
dences tMarconi et al, 2004 and Hopkins et aL 2006.) . We therefore 
present a plot of bolometric corrections against luminosity with 
these classes of objects removed in Fig. [8] These sources are also 
highlighted in the aox — Li/ (2500A) plot in Fig.|4] We now find a 
weak increase in bolometric correction with luminosity with these 
objects removed, but with far more deviation than the forms sug- 
gested in the literature. 

4.3 UV and Optical variability: Simultaneity of SED data 

The presence of AGN variability in the X-ray and UV bands and 
the l inks between variabili ty in the two bands are well documented 
(e.g. lShemmeretalJl200lh . For the AGN in this sample, SED data 
points have been taken from different publications in the literature, 
calculated from observations which are not in general contempo- 
raneous. While we might expect delays of a few days between ob- 
servations to be insignificant and in accordance with the time de- 
lays expected between X-rays and optical wavelengths, the fact that 
the observations used to construct the SEDs are often separated by 
many years implies that SED data points in different bands could be 
providing snapshots of the SED in epochs with particularly incon- 
gruous spectral shapes. This could partially account for some of the 
poorer model fits for some AGN in the sample. The lack of simul- 
taneous SED data could therefore represent a significant obstacle 
to es timating bolometric cor rections accurately. 

iBrocksoDP et al.l ilOO fh present simultaneous optical-UV-X- 
ray points for 23 AGN, using XMM Optical Monitor (XMM- 
OM) data for optical and UV points. Fourteen are common to 
our sample, allowing us to gauge the importance of simultaneity 
in determining accurate SEDs. These authors generate SED fits 
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Figure 9. Bolometric corrections generated from simultaneous SEDs (open 
blue diamonds - using XMM Optical Monitor optical-UV data and XMM 
X-ray fluxes as given in lBrocksopp et alj200 6l) and non-simultaneous SEDs 
(filled black circle s - FUSE UV points, m iscellaneous optical, XMM X-ray 
fluxes as given in lBrocksopp et al]|2006l) . Other conventions as in Fig. [3] 
The AGN in the subsample: 1 - PG 0947+396, 2 - PG 0953+414, 3 - PG 
1048+342, 4 - PG 1114+445, 5 - PG 1115+407, 6 - PG 1116+215, 7 - PG 
1216+069, 8- PG 1322+659, 9- PG 1352+183, 10 - PG 1402+261, 11 - PG 
1415+451, 12 - PG 1444+407, 13 - PG 1512+370, 14 - PG 1626+554, 15 - 
3C 273. 

to the data using a more detailed Comptonized disc model than 
the DISKPN model used in our study, but for the purposes of cal- 
culating bolometric corrections, our simp le model should be suf- 
ficient. Indeed, fitting the data points in iBrocksopp et al ] ( l2006h 
to the DISKPN and power law combination yields better fits than 
our non-contemporaneous datasets for these AGN. The bolometric 
corrections calculated from these fits are presented in Fig. |9] with 
the bolometric corrections calculated from non-contemporaneous 
datasets presented for comparison. 

We see from Fig.|9]that the bolometric corrections from simul- 
taneous SEDs trace out similar relative positions to each other as 
for non-simultaneous SEDs, but are generally higher than for non 
simultaneous SEDs. Inspection of the SEDs (Fig.llOt reveals that 
this is primarily due to the XMM-OM points (between wavelengths 
2120A < A < 5483A) missing the peak of the BBB emission in 
the UV; this allows the model to peak higher in flux than when 
constrained by the FUSE UV points. As discussed previously, the 
dip in the FUSE UV points at the peak of the BBB could be due to 
intrinsic reddening, but we do not explore this issue further here. It 
is important to note from this that we can only recording the true 
accretion luminosity if we use simultaneous SEDs; this motivates a 
further study with a larger simultaneous sample. 

4.4 X-ray Variability 

For those AGN for which multiple ASCA observations were avail- 
able, the average and standard deviation X-ray fluxes and photon 
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Figure 8. Hard X-ray Bolometric Correction against 2-10 keV luminosity, with Radio Loud, NLSl and X-ray weak AGN removed. Lower panel: binned 
bolometric connections, six AGN per bin. Other conventions as in Figure[3] The AGN in the subsample: 1 - PG 1116+215, 2 - PG 1322+659, 3 - Ton 951, 
4 - PG 0026+129, 5 - Ark 120, 6 - PG 0052+251, 7 - PG 0804+761, 8 - PG 0947+396, 9 - PG 0953+414, 10 - PG 1114+445, 11 - PG 1216+069, 12 - PG 
1307+085, 13 - PG 1415+451, 14 - PG 1444+407, 15 - PG 1626+554, 16 - PG 1229+204, 17 - Mrk 205, 18 - Mrk 79, 19 - Mrk 279, 20 - MR 2251-178, 21 
- Mrk 1383, 22 - Fairall 9, 23 - NGC 985, 24 - PG 2349-014, 25 - NGC 4395, 26 - NGC 4593, 27 - NGC 3516, 28- NGC 3227*, 29 - NGC 7469, 30 - NGC 
6814, 31 - NGC 3873, 32 - Mrk 290, 33 - PG 1352+183, 34 - ESO 141-G55, 35 - PG 1048+342, 36 - II Zw 136, 37 - Mrk 509, 38 - Mrk 506, 39 - Mrk 876. * 
See footnote for this object in Fig. [5] 
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Figure 10. Comparing a sample of AGN SEDs from simultaneous 
data (XMM X-ray and XMM-OM optical and UV points) and non- 
simu ltaneous data (FUSE U V, miscellaneous optical, XMM X-ray fluxes 
from lBrocksopp et al.l200d) . Filled cricles represent non-simultaneous data 
points, open squares represent simultaneous points. 



indicies were employed which feeds directly into the determination 
of errors on the bolometric correction, thereby taking into account 
some element of X-ray variability. In other cases, where multiple 
observations from different satellites or different publications in 
the literature were available, it is not robust to take averages due 
to potentially differing systematics. In these cases we compare the 
ratios of 2-10 keV fluxes from the different satellites/publications 
in table |2] For the 12 AGN where observations from both ASCA 
and XMM were available, we find that there is no identifiable sys- 
tematic offset between fluxes from the two satellites, with XMM 
fluxes exceeding ASCA fluxes 58 per cent of the time. 

We see that a typical value for the ratio 'maximal flux to 
minimal flux' is ~ 1.5 enabling us to estimate the effect of 
this variation on the bolometric correction. For example, con- 
sidering 3C 273 with max/min flux ratio 2.19, the bolomet- 
ric corrections generated from the differ ent X-ray datasets are 
14.7t^;^ (ASCA), 28.7 86l^-^gg (XMM --iBrocksopp et al.|[200^ 
and 20Atl 'o (XMM - Turler et alj|2006l) . This gives a ratio of 
maximum bolometric correction to minimum bolometric correc- 
tion of ~ 2. It is also worth bearing in mind that some AGN are 
known to have significantly higher variability than the subsample 
highli ghted here, a case in point being NGC 4395 JVaughan et al.l 
I2OO5I) which has been observed to rapidly vary in X-rays by fac- 
tors of 10. There is therefore an obvious tendency for variability to 
increase the spread in bolometric corrections, and this adds to the 
uncertainties in their determination. 



4.5 Effects of Variations in SMBH mass estimates 

The normalisation of the disk component for each of the AGN in 
the sample is constrained using the mass and luminosity distance 
to the source. The bolometric corrections presented so far do not 
reflect uncertainties in SMBH mass estimates; accounting for such 
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Figure 11. Effect of varying the SMBH mass on the SEDs generated. Red 
dot-dashed lines show the model fit with mass increased by a factor of 3, 
blue dashed Unes show the fit with the mass decreased by the same factor; 
other conventions as in Fig.E] For NGC 3873, increasins the mass decreases 
K2- lOkcV by 8 per cent and decreasing the mass increases it by 40 per cent. 
For Mrk 478, increasing the mass reduces K2-iOkcV by 35 per cent and 
decreasing the mass increases it by 107 per cent. These two objects show 
particularly large changes; the average variations were significantly lower 



uncertainties is difficult due to the diverse methods used to obtain 
these estimates. 

The well established technique of reverberation mapping, for 
example, allows the mass to be estimated by tracking variations in 
the continuum flux from the accretion disc and comparing them 
with variations in emission lines from the broad line region (BLR) 
surrounding it. The delay between variations allows the size of the 
BLR to be determined, which can be used to get a virial mass es- 
timate. This technique forms the basis for the mass determinations 
for many of the AGN in the sample. Other techniques derived from 
this exploit empirical correlations between, for exa mple, optical lu- 
minos ity at 5IOOA and the BLR radius. The study of lPeterson et al.l 
( I2OO4I) suggests that there is a systematic uncertainty of a factor 
of ~ 3 in mass estimates derived using the reverberation map- 
ping method. There are also potential systematics in the estimation 
of NLSl masses, possibly due to the flattening and inclination to 
the line of sight of the broad-line region; this could cause NLSl 
BH mas ses to be underesti mated under the reverberation mapping 
method jCollin et al.ll2006l) . 

We attempt to assess the importance of such variations by 
varying the SMBH masses up and down by a factor of 3, since this 
represents a typical systematic uncertainty in BH mass estimation 
methods. We exclude MR 2251-178 since the m ass estimate used is 
already an upper limit ^Morales & Fabiatj|2002l) . A sample of SED 
fits generated with the varying BH masses are shown in Fig.ll II 

We find that if the mass is increased by a factor of 3, the bolo- 
metric correction decreases by 14 ± 13 per cent on average, and 
decreasing the mass by the same factor increases the bolometric 
correction by 8 ± 15 per cent on average. If the fits produced from 
varying the mass produced reduced chi-squared values greater than 
10 times higher than their previous values, we discarded these ob- 
jects from the average since these are likely to represent particularly 
poor fits to the data giving unrealistic SED shapes. The averages 
presented are of a similar order to the random error in bolometric 
corrections due to other sources, and so systematic uncertainties of 
this order (a factor 3) in the mass are significant, but not more so 
than uncertainties from other sources. 
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AGN 


Available X-ray observations [ 1 ] 


Observation with maximal 
2 - lOkcV flux 


Maximal 2 

(10-12 


- lOkeVflux 

erg/s/cm^) 


Maximal flux/Minimal flux 


3C 273 


ASCA (T), 

AiVliVi yij ), 

XMM (Turler et al. 2006') 


ASCA (T) 




98.91 


2.19 


IZw 1 


ASCA (T), 

AiVliVi yr) 


XMM (P) 




4.83 


2.86 


Mrk 478 


ASCA (T), 

AiVliVi \r) 


XMM (?) 




3.18 


1.59 


PG 0947+396 


XMM (B), 

A.iVliVi yr) 


XMM (P) 




1.73 


1.18 


PG 0953+414 


ASCA (T), 
XMM (B), 

AiVliVl \r) 


XMM (?) 




2.80 


1.26 


PG 1048+342 


XMM (B), 

AlVliVl yr) 


XMM (B) 




1.37 


1.15 


PG 1 100+772 


ASCA (T), 

AiVliVl yr) 


XMM (?) 




3.07 


1.67 


PG 1 1 14+445 


ASCA (T), 
XMM (B), 

YAylA-T /P\ 
AiVliVl yr) 


XMM (B) 




2.43 


1.09 


PG 1115+407 


XMM (B), 
XMM (P) 


XMM (?) 




1.21 


1.10 


r\j III 0+Z i J 


yl ), 
XMM (B) 


A A /T'l 

yi) 




J. i\j 






yl), 
XMM (B), 

AiVliVl yr) 


A C/^ A /"T\ 

AoCA yL) 




1 '\A 


on 


PG 1322+659 


ASCA (T), 

AiVliVl 1^13 ), 

XMM (P) 


XMM (?) 




1.12 


1.66 




YA/fA-T iTi\ 

AiVliVl yt) ), 
XMM (P) 


YA/fA/T ^'P^ 
AiVliVl yr) 




1 SQ 
i.oo 






ACr'A (T\ 
/\ . 3 V r\ yL), 

XMM (B) 


/\ . 3 V r\ yi) 




1.7 J 


l.f J 


Pn 1 Al S-uA'^ 1 


AiVliVl yij ), 
XMM (P) 


AiVliVl yr) 




1 C\f\ 
i.UO 


i.U/ 


PG 1444+407 


ASCA (T), 
XMM (B) 


ASCA (T) 




0.59 


1.41 


PG 1512+370 


XMM (B), 
XMM (P) 


XMM (?) 




1.87 


1.11 


PG 1626+554 


XMM (B), 
XMM (P) 


XMM (B) 




3.12 


1.03 


TON 951 


ASCA (T), 
XMM (P) 


XMM (?) 




3.93 


2.33 



Table 2. Hard X-ray variability for AGN for which multiple satellite observations were availa ble. [1] ASCA (T) denotes ASCA observations from the 
Tarta rus database, XMM (B) denotes XMM observations as presented in lBrocksopp et all i2006ll and XMM (?) denotes XMM observations as presented in 
IPiconcelli et al. 12005). 



5 PARALLELS WITH GALACTIC BLACK HOLE 
ACCRETION 

Our results may provide another perspective on the question of 
whether AGN accretion processes can be understood as a scaled-up 
versio n of Galactic Black Hole (GBH) accretion. McHardv et al. 
(l2006h discuss the importance of X-ray emission in this regard, 
since it emerges from very close to the black hole. They specif- 
ically discuss the topic of a universal relation between black hole 
mass and 'break timescale' (timescale at which power spectral den- 
sity steepens) across B H masses ranging from Stellar to Super- 
massive. The review o f iRemillard & McClintockl (l2006h discusses 
the three state picture for GBH accretion, identifying the 'hard' 
state (formerly known as the low/hard state), the thermal state (for- 



merly the 'high/soft' state) and the steep power-law (SPL) state 
(formerly the 'very high' state); along with the preponderance for 
quasi-periodic oscillations in the thermal/SPL states. The Edding- 
ton ratio Lboi/LEdd is a useful separator of these states, and there 
is a progression from low to high Eddington ratios as these sources 
undergo transitions between the hard, thermal and SPL states. The 
latter is often observed close to the Eddington limit. In this sec- 
tion, we discuss our results in the context of the emergent unified 
pictures of accretion on all mass scales. 



5.1 Bolometric Correction as a function of Eddington Ratio 

We present the bolometric corrections for AGN plotted against Ed- 
dington ratio in Fig. [12] top panel, with the outlier PGlOl 1-040 
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Figure 13. SEDs (normalised at 1 eV) for AGN with Eddington ratios in 
the range 0.0012 < Lboi/LEdd < 0.032 (black solid curves) and 0.27 < 
Lbol/LEdd < 2.7 (blue dashed curves). The thicker curves show the mean 
SEDs calculated from the two sets of SEDs. 

(X-ray weak) removed. The Eddington ratio can be understood as 
the luminosity scaled by the mass of the central black hole, and 
presented in this form there appears to be a distinct step change 
in bolometric correction at an Eddington ratio of ~0. 1. We con- 
firm previous results which highlig ht the tendency f or NLSl nu- 
clei to have high Eddington ratios (.Pounds et al.ll99S ). If accretion 
states in AGN are qualitatively similar to those in GBHs, the X-ray 
bolometric corrections could be a useful diagnostic of the accre- 
tion state in a black hole and the step-up in bolometric corrections 
above ~ 0.1 Eddington could be indicative of AGN undergoing a 
tran sition between such accretion states. 

IWang et al] | |2004|) present the variation of L2-iokcv/Lboi 
(i.e. l/K2-iokev) against Eddington ratio for a sample of 
AGN. Their their results agree qualitatively with ours; how- 
ever, their bolometric luminosities are not always calculated 
from the SED continuum, as some of them are estimated us- 
ing Lboi ~ 9ALa(5100A) jKaspi et all 120001) . For those objects 
where bolometric luminosities are integrated from the multiwave- 
length continuum, they use values from ,Woo & Urrvl j2002t) and 
iGrupe et al.l ( l2004l) . who do not consistently integrate over the same 
ranges (the former integrate from infrared to X-rays, whereas the 
latter integrate from optical to X-rays as we do here). It is there- 
fore possible that our approach of constructing SEDs afresh and 
integrating over the same energy range in each case may provide a 
window on some of the intricacies in the variation of the bolometric 
correction with Eddington ratio (such as the step change perceived 
here). 



5.1.1 Accretion States 

Previous work has allowed the characterisation of distinct spec- 
tral features in the different accretion states of GBH binaries. At 
low Eddington ratios, the low/hard state is characterised by a hard 
(r < 2) power law, often dominating the black body disk spec- 
trum. In the high/soft state, the power law spectrum is much softer 
(r > 2), but often the power law again dominates over the disk 



component. In intermediate states such as the thermal state, the 
black body is comparatively more prominent. There is a progres- 
sion for increasing overall flux in the various components going 
from low to high states. To contrast these characteristics with the 
AGN SED at high and low Eddington ratios, we reproduce the 
SEDs normalised at O.I eV, for eight AGN in each of the bins 
for Eddington ratios 0.0012 < Lboi/LEdd < 0.032 (bin I) and 
0.27 < Lboi/IjEdd < 2.7 (bin 2) in Fig. [T3](we do not use the 
highest Eddington ratio bin since uncertainties surrounding mass 
determinations could be artificially producing extremely super- 
Eddington accretion rates). 

The average SEDs for the high and low Eddington AGN al- 
low us to draw some parallels with high and low accretion states 
in GBH binaries. The lower Eddington AGN have a harder spec- 
tral shape in X-rays, as expected for the low/hard state, and the 
higher Eddington AGN have softer X-ray spectra as expected of 
the high state. The dependence of the X-ray photon index on Ed- 
dington r atio in AGN has been commented on by other authors 
including lshemmer et al. IJiooi), who find a correlation between F 
and Lboi/IjEdd. Their findings are broadly echoed by our average 
SEDs. In this study, we are also able to see how optical-UV emis- 
sion varies with Eddington ratio. We find that the shift in the peak 
energy of the SED is roughly that expected from the dependence 
of Eddington ratio on temperature; For a disk the black body lu- 
minosity scales as L cx T^M^, and since LEdd oc M, The ratio 
L/Lsdd will scale as T^M. The two bins used span a factor of 
~ 140 in Lboi/IjEdd; the expected range in temperatures would 
therefore be a factor of ~ 3.4, assuming a similar distribution of 
masses for the AGN in the two bins. The ratio of temperatures ob- 
served between the two bins is indeed ~ 3.4. The assumption of 
similar mass distributions does appear to be valid here: the average 
values of log(MBH/M0) in the low and high Eddington ratio bins 
are 7.34 ± 0.38 and 7.76 ± 0.19 respectively. The similar mass 
distributions seen in the two bins may imply that the disk structure 
is probably similar in the 'hard' and 'soft' states, and is broadly 
consistent with the hypothesis that the disk extends inwards close 
to the black hole without truncation. 

5.1.2 Coronal Fraction 

We can use the 2-10 keV bolometric correction to calculate the 
total coronal fraction for our two fiducial high and low Eddington 
AGN SEDs, by extrapolating the coronal emission over the whole 
X-ray energy range. We employ the mean photon indices for high 
and low Eddington ratios to calculate the factor difference between 
the 0.1-250 keV flux and the 2-10 keV flux in the two states. We 
can then estimate the power-law fraction, or 'coronal fraction', us- 
ing typical bolometric corrections of ^ 20 for low Eddington and 
~ 50 for high Eddington. The observed coronal fractions emerge 
as ^ 35 per cent for low Eddington ratio and ~ 10 per cent for high 
Eddington ratio, representing a fall of almost a factor of 4 between 
low and high Eddington ratio regimes. 

However, the presence of X-ray reflection introduces some 
subtle modifications to this picture We have so far assumed that 
the X-ray emission represents the total coronal emission Lc and 
the UV bump (LD)is solely from the accretion disk. If reflection 
is taken into account, the X-rays we see are in fact only part of 
the total coronal emission (assumed to be 1/2 for simplicity here, 
corresponding to a reflection fraction of unity), and the rest is radi- 
ated onto the disk. The UV bump then contains both emission from 
the thermal accretion disk and some fraction / of reprocessed X- 
rays from the corona (Lroproc = 0.5fLc). The true coronal fraction 
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Figure 12. Upper panel: 2-10 keV bolometric correction against Eddington ratio for the AGN in the sample. Blue empty squares represent NLSl nuclei, 
black empty circles represent radio loud AGN and dark red points with enor bars are binned results. Lower panel: bolometric corrections for the GBH GX 
339-4 (red triangles) superimposed on binned AGN results. Individual AGN data points removed for clarity. Grey shaded region shows the interquartile range 
calculated from Monte-Carlo realisations of the model fits, and the grey line shows the median. The possible problems due to the inclusion of Lboi on both 
axes are addressed in section 15.1.41 
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Figure 14. Plot of aox against L^oi/LEdd for the sample. Key as in Fig. 


Ftruo = Lc / (Lc + Ld) will then differ from the observed fraction 
Fobs = 0.5Lc /(0.5Lc + Ld + 0.5£Lc). Using a simple pexrav 
reflection model in XSPEC with reflection fration as above (canon- 
ical abundances assumed), we obtain absorption fractions f ~ 2/3 
and ~ 5/6 for low and high Eddington photon indices respectively. 
This results in true coronal fractions of 45 per cent for low Edding- 
ton and 11 per cent for high Eddington. 

5.1.3 Implications for Surveys 

iKoUmeier et al. I (12006!) find a dearth of AGN below Lboi/LEdd ~ 
0.1 in a joint X-ray - IR - optical study of the Bootes field. If 
objects are indeed missing then they would have masses of about 
10^ Mq, and so are not well represented by those in our study. If 
there is an abrupt change in bolometric correction of these objects 
at an Eddington ratio of 0.1 such as found here, then this may in- 
fluence their selection of objects. Our results do not obviously ac- 
count for their finding in a straightforward way, although we note 
that objects in the high and low Eddington ratio bins shown in Fig. 
[T3]have very different ionizing photon fluxes (here taken to be the 
flux between 0.036 and 0.1 keV). Using the SEDs from Fig. 1131 
we calculate that the higher Lboi /Lsdd AGN have ~ 40 times the 
ionizing photon flux of the lower Lboi/LEdd AGN, despite hav- 
ing similar optical to X-ray ratios (aox). This could affect the ob- 
served line luminosities and hence also any selection criteria based 
on them. We find further confirmation of this in Fig. [141 where we 
present aox against Lboi/IjEdd. The relatively constant value of 
qqx seen across all Eddington ratios in the sample reiterates the 
point that the optical-X-ray ratios in AGN do not provide us with 
sufficient information on the important variation in the UV which 
is needed for calculating accurate bolometric corrections. We see 
from Figs.[T3]and[T2]that this variation is significant. 

We can briefly comment on the traditional estimation tech- 
nique of using 9AL a(5100A) to e stimate the bolometric luminos- 
ity, taken from Ka spi et alj feOOOl ). The optical bolometric correc- 
tion Lboi/ALA(5100A) is plotted against Lboi/LEdd in Fig. [TS] 
We notice that for low Eddington AGN, the bolometric correction 
of 9 (with a scatter comparable to the factor ~ 2 identified by 
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Figure 15. Optical bolometric correction Li3oi/AL;j(5100A) against Ed- 
dington ratio Lboi/IjEdd. Key as in Fig. ll2l (the red diamonds represent 
X-ray weak objects as i n previous figures ). The dotted line shows the bolo- 
metric correction of 9 f Kaspi et al.l l2000) which is often used to estimate 
bolometric luminosities from monochromatic optical luminosities. 

[Richards et ^|2006|) is relatively accurate, but deviates from this 
at high Eddington ratios. The optical bolometric correction shows 
a similar trend to the X-ray bolometric correction, increasing with 
Eddington ratio. Deviations from the traditional estimation tech- 
nique may therefore be particularly important for NLS Is, if indeed 
their super-Eddington accretion rates are real. However, since the 
optical points may be prone to some uncertainties such as contam- 
ination from host galaxy starlight etc., we caution against extrap- 
olating too much from this plot, apart from the general trend for 
increasing bolometric correction with Eddington ratio. 



5.1.4 Consideration of Systematics 

It is important to note that where bolometric correction is plotted 
against Eddington ratio, there is the potential to introduce false cor- 
relations since Lboi is present on both axes. However, if we gener- 
ate a simple scatter plot using randomly generated input parameters 
taking values in the ranges available in our data, we find the slope 
of any implied or 'false' trend is far steeper than the trend seen 
here, leading us to believe that the trends we see are real. However 
the lack of objects in the lower-right regions of the bolometric cor- 
rection against Lboi/LEdd plots could be in part due to this. There 
is also the issue of the preponderance for high mass objects to have 
high X-ray luminosities. Neither of these systematics would repro- 
duce the 'step' change in bolometric correction seen here, and so 
we propose that these systematics are not immediately relevant. 

Finally, we note the possibility that the selection bias inherent 
in the FUSE sample could be artificially producing the step-change 
seen here at ~ 0.1 Eddington. The AGN represented in the FUSE 
sample are relatively UV bright and also low in redshift. The large 
number of AGN clustering at Eddington ratios of ~ 0.1 could, for 
example, reflect a tendency for local AGN to have this accretion 
rate, explaining in part the lack of significant numbers of low Ed- 
dington AGN. The UV brightness of the sample would reinforce 
this, implying a bias away from 'low state' AGN, which are ex- 



© 0000 RAS, MNRAS 000, 000^00 



Bolometric Corrections for AGN 1 5 



pected to be UV-faint (Figure [T3t. The domination of NLSl nuclei 
at iiigh Eddington may also not be representative of the popula- 
tions making up the XRB. However confirmation or refutation of 
this will have to wait for a more extensive study including AGN at 
higher redshifts. 

5.2 Comparison with GX 339-4 

In order to compare the bolometric corrections between AGN and 
GBHs, we need to bear in mind how different components of the 
emission scale with mass; the scaling of the disk temperature with 
mass causes the 'BBB' in the UV arising from the disk emission 
in AGN to be shifted to X-ray energies in GBHs. However, the 
power-law component arising from the coronal emission is scale- 
independent. So to estimate the total accretion flux in GBHs we 
evaluate the flux from the lowest energy of the disk bump until 
^ lOOkeV, and then to evaluate the 2-10 keV power-law flux we 
neglect the disk component to obtain only the coronal flux in 2- 

10 keV. We use fits to RXTE, Integral and Beppo S AX data on the 
GBH GX 339-4, provided by 'Miller et alj j2006h . ICorongiu et al.l 
( l2003h (ha rd state),[B elloni et al. (2006) (soft intermediate/thermal 
state) and lMiller etal.1 <2004l) (SPL state) to calculate bolometric 
corrections. The 2-10 keV bolometric corrections against Edding- 
ton ratio are presented in the lower panel of Fig. [12] 

It is interesting to note that the bolometric corrections do in- 
crease as expected as the source moves between the hard, thermal 
and SPL states (increasing Eddington ratio). The transitional region 
between the hard and SPL states is at around Lboi/LEdd ~ 0.1, 
allowing us to tentatively identify the high and low bolometric 
correction regions seen in AGN as being from two relatively dis- 
tinct accretion states. It is known from the literature, however, 
that a number of GBH binaries undergo state transitions at signifi- 
cantly lower Eddin gton ratios, of the order ~ 2 per cent Eddington 
( lMaccaron3 l2003h . They discuss the variation in transition lumi- 
nosity in going from hard-to-soft and soft-to-hard states, and the 
hysteresis implied by it. It is possible that we may be witnessing 
similar hysteresis in our AGN sample by virtue of the lack of a 
very sharp, 'thin' transition between low and high bolometric cor- 
rections around Lbo i/LEdd ~ O-l- 

The work of iKording et all (2006^ discusses the use of 
'Hardness-Intensity' or 'Disk Fraction Luminosity' diagrams 
(HIDs or DFLDs) for characterising the accretion states in both 
AGN and X-ray binaries. To construct the AGN HID, They define 
the 'power-law luminosity' LpL as 3 x L2-iokoV and calculated the 
disk luminosity Ld by applying a bolometric correction from the 
B-band. The Hardness Ratio (HR) is defined as Lpl/(Ld + Lpl). 
The HID consists of Ld + Lpl plotted against the HR. We at- 
tempt to construct a HID for our sample, but here we make use 
of the UV SED fit to work out the disk luminosity directly, defining 
Ld = Lo. 001-0. ikcv- The resulting HID is presented in Figure [76l 

We note firstly that there are no objects in our sample with 
0.8 < HR < 1.0. The lack of a significant 'hard' population 
means that we can only make a limited comparison with their work 
using our sample. However, we note that the high Eddington NLSl 
nuclei are clustered in the region expected for soft-state objects, and 
the r adio loud objects also lie where expected as given by Figure 

11 of lKording et"al] ( l2006h . 

GBH binaries show striking changes in rad io flux between 
hard and soft states (e.g. iFender & Bellonill2004|) . The radio be- 
haviour of AGN is less clear. There have been a number of recent 
studies of the radi o loudness of AGN versus E ddington ratio (e.g. 
ISikora et al.ll2006l . IPanessa et alj|200l iMaoj [2007.) which show 




0.2 0.4 0.6 0.1 

Hardness Ratio = Lpl / (Ld + Lpl) 

Figure 16. HID/DFLD for the AGN in our sample. Key as in Fig.fTD 



that higher radio loudness correlates with with lower Eddington 
fractions, but no sharp switch is se en. Using 5 GHz ra dio lumi- 
nosities from NED and the study of iMerloni et al.l ( l2003h . we find 
that in our small sub-sample of 16 AGN in bins 1 and 2, there is 
no identifiable tendency for these low Eddington AGN to be radio 
loud (using the ratio of radio to X-ray luminosity as a measure of 
radio loudness). A recent study bv lXue & Cuil j2007h discusses the 
possibility that even the correlations seen in GBH binaries may be 
influenced by significant systematics. The radio - X-ray connection 
in BH accretion therefore still requires clarification. 



6 SUMMARY AND CONCLUSIONS 

We have demonstrated that the hard X-ray bolometric correction 
for AGN may vary significantly from object to object, especially 
when considering classes of AGN such as Narrow Line Seyfert 1 
nuclei and Radio Loud AGN. Removing these classes of objects 
yields only a weak trend of increasing bolometric correction with 
luminosity (Fig.[8l(. The uncertainties introduced by the lack of si- 
multaneous SEDs, X-ray and systematics in SMBH mass estima- 
tion methods are also significant. We a rgue that the d ependences 
of K2-iokeV on L2- iokcV proposed by iMarconi et al . (2004) and 
Hopkins et al.l j2006l) do not reflect the spread observed in the real 
AGN population because they primarily take into account the SED 
shape via the aox — L^ (2500A) correlation, missing the intrinsic 
spectral differences attributable to mass and luminosity seen at the 
peak of the BBB at around ~ lOOOA. We show that differences 
are important and are the main cause of the spread in bolometric 
corrections. 

At first sight this could be accounted for by invoking intrinsic 
reddening in UV-optical part of the SED, but the available infor- 
mation on reddening seems to suggest that this will increase bolo- 
metric corrections by at factors of ~ 2 at most, without necessar- 
ily aligning the bolometric corrections with the suggested forms. 
There may also be a luminosity dependence of intrinsic redden- 
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ing which could go against the trend for increasing K2_iokeV with 
luminosity. An understanding of intrinsic extinction therefore re- 
mains of paramount importance and further study into its effects 
in all types of AGN would be extremely valuable. However, we 
establish that the effect of any of reddening corrections will be to 
increase the X-ray bolometric corrections. In order to satisfy the 
Soltan argument, this would require an increase in the mass-to- 
radiation conversion efficiency, 77. Higher efficiences could poten- 
tially imply more rapidly spinning black holes and de-emphasize 
the relative importance of radiatively inefficient modes of accretion 
in the X-ray Background. 

There are some potential limitations to this study due to the 
SED construction method. These findings rely on the fitting of 
multiband data to a simple black-body disk model and power law. 
In order to get a more robust and physically realistic SED and to 
improve so me of the poor fits, m ore advanced models such as those 
discussed by Blaes et al. I( l200ll) for 3C 273 could be employed. The 
use of only four to eight data points may not sufficiently constrain 
the SED fit to give an accurate bolometric correction, and the cur- 
rent method may miss out potentially important contributions to the 
accretion luminosity, such as the soft X-ray excess. However the 
good agreement between aox values from our SEDs and the sam- 
ple of Strateva et al. ( 2005) would imply a sensible reconstruction 
of the continuum SED shape for most of the AGN in this sample. 
Additionally, there is still much uncertainty as to how well the cur- 
rent disk models fit to the UV disk emission. With the absence of 
any clear picture as to which is the 'correct' model to be used for 
the AGN disk emission, the DISKPN model used here provides a 
sensible first order approach to estimating the UV emission. It is 
also clear that the extreme bolometric corrections (both high and 
low) seen in our sample are independent of such effects and so the 
(hitherto unseen) degree of spread in K2-iokcV seen in our results 
can be confirmed as real. 

While the relation between bolometric corrections and lumi- 
nosities in the 2-10 keV band seems subject to significant uncer- 
tainties, it appears that Eddington ratio is a better discriminator of 
populations of AGN with high and low bolometric corrections. We 
find a distinct tendency for bolometric corrections of 15 — 25 below 
a threshold of around Lboi/LEdd ~ 0.1, with values of 40 — 70 at 
higher Eddington ratios. Making allowances for reflection, we find 
that this corresponds to a fall in the coronal fraction of a factor of 
~ 4 when transitioning from low to high Eddington regimes. The 
transitional region at Lboi/LEdd ~ 0.1 is of particular interest be- 
cause of the potential parallels to be drawn with GBH accretion 
states. Notably, bolometric corrections for the GBH candidate GX 
339-4 calculated from observations in the literature also display a 
transition at around Lboi/LBdd ~ 0.1, marking the transitional re- 
gion between hard and soft states. However, effects of hysteresis 
in state transitions on both SMBH and GBH scales could blur this 
simple picture. Additionally, the aforementioned uncertainties due 
to intrinsic reddening, UV/optical and X-ray variability and sys- 
tematic errors in mass estimates may again introduce uncertainties, 
and selection biases inherent in our sample must also be consid- 
ered. 

We note in passing that nearly one half 45 per cent) of 
the power of low Eddington ratio AGN appears to emerge in the 
corona, indicating that much of the accretion flow is extracted by 
magnetic fields and dissipated outside the disk (taking X-ray re- 
flection into account). For high Eddington ratio AGN, this fraction 
reduces to approximately one tenth. 

The significant differences between bolometric corrections 
from non-simultaneous and simultaneous SEDs implies that to gain 



an accurate understanding of the AGN SED, a much larger cata- 
logue of simultaneous SEDs akin to the sample of lBrocksopp et all 
( 2006) is a pressing need. 

We finally note that we are presented with an opportunity to 
refine X-ray background modelling techniques further than has cur- 
rently been done. The pioneering studies of previo us authors started 
off wi th an assumption of constant K2-iokcV (e.g.lElvis et al.'200i 
iFabia n & Iwasawa 1999), followed by the studies of Fabian (2003) 
and Marconi et ah) ( 120041) who pointed out the importance of tak- 
ing into account a luminosity dependence. In particular, the lat- 
ter study illustrates how using a physically realistic dependence of 
«:2-iokcV on luminosity yields a SMBH mass density from the X- 
ray background which agrees with that determined by matching the 
local and relic SMBH mass functions. This work suggests that the 
more fundamental dependence may be on Eddington ratio, and a re- 
formulating the modelling process to incorporate this dependence 
may be useful for future studies. 
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AGN name 


Redshift 


Optical data source 


UV data source 


X-ray data source 


Mass data source 


NGC 4395 


0.001064 


(None found) 


[1] 


ASCA 


[2] 


NGC 3227 


0.003859 


HST [3] 


HST-STIS [3] 


ASCA 


[4] 


NGC 6814 


0.006 


Hale telescope ([5], [6]) 


lUE [5] 


ASCA 


[7] 


NGC 3516 


0.008836 


HST [8] 


HST-STIS [8] 


ASCA 


[4] 


NGC 4593 


0.009 


(None found) 


lUE [9] 


ASCA 


[4] 


NGC 3873 


0.01 


[10] 


FUSE 


ASCA 


[11] 


NGC 7469 


0.016 


Mt. Wilson telescope [12] 


FUSE 


ASCA 


[11] 


Mrk 79 


0.022 


(None used - [10] available) 


FUSE 


XMM [13] 


[11] 


Ark 564 


0.025 


(None used - [10] available) 


FUSE 


ASCA 


[14] 


Mrk 335 


0.026 


[15] 


FUSE 


ASCA 


[11] 


Mrk 279 


0.03 


[10] 


FUSE 


ASCA 


[16] 


Mrk 290 


0.03 


HSTandKPNO[17] 


FUSE 


ASCA 


[15] 


Ark 120 


0.032 


[10] 


FUSE 


ASCA 


[11] 


Mrk 509 


0.034 


HST andKPNO[17] 


FUSE 


ASCA 


[11] 


ES0141-G55 


0.036 


[10] 


FUSE 


ASCA 


[11] 


KUG 1031+398 


0.042 


(None found) 


FUSE 


ASCA 


[18] 


NGC 985 


0.042 


[10] 


FUSE 


ASCA 


[4] 


Mrk 506 


0.043 


HST andKPNO[17] 


FUSE 


ASCA 


[19] 


Fail-all 9 


0.047 


(None used - [10] available) 


FUSE 


ASCA 


[11] 


3C 382 


0.058 


(None used - [20] available) 


FUSE 


ASCA 


[21] 


PG 1011-040 


0.058 


[15] 


FUSE 


ASCA 


[15] 


IZw 1 


0.061 


[15] 


FUSE 


ASCA 


[15] 


TON SI 80 


0.062 


UK Schmidt telescope [22] 


FUSE 


ASCA 


[4] 


II Zw 136 


0.063 


[15] 


FUSE 


Ginga [23] 


[15] 


PG 1229+204 


0.063 


[15] 


FUSE 


ROSAT [24] 


[15] 


TON 951 


0.064 


HSTandKPNO[17] 


FUSE 


XMM [25] 


[15] 


MR 2251-178 


0.066 


[26] 


FUSE 


Chandra [27] 


[28] 


Mrk 205 


0.071 


(None used - [26] available) 


FUSE 


ASCA 


[29] 


Mrk 478 


0.079 


[15] 


FUSE 


ASCA 


[15] 


PG 1211+143 


0.081 


[15] 


FUSE 


ASCA 


[11] 


Mrk 1383 


0.086 


[15] 


FUSE 


ASCA 


[15] 


PG 0804+761 


0.1 


[15] 


FUSE 


ASCA 


[11] 


PG 1415+451 


0.114 


[15] 


FUSE 


XMM [30] 


[15] 


Mrk 876 


0.129 


[15] 


FUSE 


XMM [30] 


[15] 


PG 1626+554 


0.133 


[15] 


FUSE 


XMM [30] 


[15] 


PG 0026+129 


0.142 


[15] 


FUSE 


ASCA 


[11] 


PG 1 1 14+445 


0.144 


[15] 


FUSE 


XMM [30] 


[15] 


PG 1352+183 


0.152 


[15] 


FUSE 


XMM [30] 


[15] 


PG 1 1 15+407 


0.154 


[15] 


FUSE 


XMM [30] 


[15] 


PG 0052+251 


0.155 


HSTandKPNO [17] 


FUSE 


ASCA [31] 


[15] 


PG 1307+085 


0.155 


[15] 


FUSE 


XMM [30] 


[11] 


3C 273 


0.158 


HST andKPNO[17] 


FUSE 


ASCA 


[11] 


PG 1402+261 


0.164 


Univ. of Hawaii 2.2m [32] 


FUSE 


ASCA 


[33] 


PG 1048+342 


0.167 


[15] 


FUSE 


XMM [30] 


[15] 


PG 1322+659 


0.168 


HST andKPNO[17] 


FUSE 


XMM [34] 


[15] 


PG 2349-014 


0.174 


HST andKPNO[17] 


FUSE 


ASCA 


[33] 


PG 1116+215 


0.176 


[15] 


FUSE 


ASCA 


[33] 


4C +34.47 


0.206 


HSTandKPNO [17] 


FUSE 


XMM [35] 


[36] 


PG 0947+396 


0.206 


HST andKPNO[17] 


FUSE 


XMM [30] 


[15] 


PG 0953+414 


0.243 


HST andKPNO[17] 


FUSE 


ASCA 


[15] 


PG 1444+407 


0.267 


[15] 


FUSE 


ASCA 


[33] 


PG 1 100+772 


0.311 


HSTandKPNO [17] 


FUSE 


ASCA 


[15] 


PG 1216+069 


0.331 


[15] 


FUSE 


XMM [30] 


[15] 


PG 1512+370 


0.371 


[15] 


FUSE 


XMM [30] 


[15] 



Table Al. Data sources for Optical, UV, X-ray and SMBH mass estimates for AGN in the sample. FUS E observations are all t aken f rom I Scott et al ] j2004) 
and ASCA observations are from t he TARTARUS databa se, unl ess othervyise stated. Other references: [11 - [ C renshaw et al.l j2004^. [2] -'P etersoii et^i2005l) 
[31 -ICrenshaw et al.l feOOlh, [41 -lO'Neill et al.l j2005l). [51 -ICzernv et al.l il995l). [61 -Ide Bruvn & Sargenll il978l) [7] ^ aor 12001) . [81 - [Edelson et al] 



bOOdl. [91 -ISantos-Lleo et alJ h994. [1 01 -'de Vaucouleurs et al.' h99t). [1 11 -'Kaspi et al." (200&). [121 -'Anderson' (1970). [13] - Gal lo et all 
- iRomano et a l j ^20041). [ 15_l jBaskin &~ Laor (2005). [16] - Peterson et al. (2004), [17] - Shang et al. (2005), [18] - Bian & Zhao (2004) , [19] - 



.200; 



[14] 



1983 



[24] -ICeballos & Barcons ( 1996 
[291 -Bonning et al. 12005), [30] ■ 



[20] - Smith & Heckman (1989), [21] - Wu & Liu (2004), [22] - Maddox et al. ( 1990), [23] - Lawson & Turner ( 19971 

[25] - Brinkmann et al. (2003), [26] - McAlaiy et al. ( 1983), [27] - Gibson et al. ( 2005) , [28] - Morales & Fabian (2 00^ . , . , , - - . . 

iPiconcelU e t al. (2005), [31] - Ued a et aTJi200l . [32] - ISurace et al.H200ll) , [33] - iMcLure & Dunlod j200lh , [34] - iBrocksopp et alj 120061) , [35] - IPage et alj 
i2004 . [361 - ,Wang et al., ( ,2004.) 
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